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ABSTRACT 

Wang et al. (2006) estimated an average radiative efficiency of 30%-35% 
for quasars at moderate redshift. We find that their method is not indepen- 
dent of quasar lifetimes and thus that quasars do not necessarily have such 
high efficiencies. Nonetheless, it is possible to derive interrelated constraints 
on quasar lifetimes, Eddington ratios, and radiative efficiencies of supermassive 
black holes. We derive such constraints using a statistically complete sample of 
quasars with black hole mass estimates from broad Mg II, made both with and 
without the radiation pressure correction of Marconi et al. (2008). We conclude 
that for quasars with L/LEdd ^ 0.02, lifetimes can range from 140 to 750 Myr 
for Schwarzschild black holes. Coupled with observed black hole masses, quasar 
lifetimes of <140 Myr would imply that radiatively inefficient accretion or BH 
mergers must be important in the accretion history of quasars. Given reasonable 
assumptions about the quasar population, if the average quasar lifetime is < 1 
Gyr, and if radiatively inefficient accretion is unimportant, then not many BHs 
with Eddington ratio < 0.2 can be rapidly spinning. 

Subject headings: black hole physics — quasars: general 



1. Introduction 



A quasar is powered by matter accreting onto a supermassive black hole (e.g., Rees 
1984). Gas orbiting in the innermost stable circular orbit (ISCO) may be perturbed and fall 
into the black hole (BH), adding mass and angular momentum to it. To have reached the 
ISCO through a thin accretion disk, gas must have radiated away a fractional binding energy 
per unit rest mass ~ 1 — (1 — 2GM/3c 2 RiscoY^ 2 which is the system's radiative efficiency 
r) (Bardeen 1970). Since Risco decreases from QGM/c 2 for Schwarzschild (nonrotating) 
BHs to GM/c 2 for co-aligned accretion onto Extreme-Kerr (maximally rotating) BHs, more 
energy is produced by (co-aligned) thin disk accretion of a given mass onto a rotating BH 
than onto a non- rotating BH (Carter 1968). 
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If a quasar accretes gas with a fixed angular momentum vector at a fixed mass accretion 
rate M, the BH spin will increase to a theoretical maximum and the luminosity L of the 
quasar will increase along with it, since L oc rjMc 2 . In many realistic accretion models, 
the spin of a supermassive BH increases rapidly and then fluctuates around a maximum 
value, although spin-down is also possible (Volonteri 2007). Several such models, including 
numerical ones by Di Matteo et al. (2005), suggest that the accretion process stops when 
a BH becomes massive enough to support a kinetic and/or radiative luminosity capable of 
blowing away the gas fuelling it (Silk & Rees 1998; Fabian 1999; King 2003). However, the 
radiative luminosity for a given mass accretion rate depends upon the radiative efficiency 
A rapidly rotating BH can help shut down the accretion process earlier than a nonrotating 
BH could. Knowledge of supermassive BH spins is therefore useful in constraining models 
of quasar development. 

Observationally, the quasi-periodic variability detected in Sgr A* may be evidence of 
rapid spin of the Galactic BH (e.g., Genzel et al. 2003). Evidence for rotating supermassive 
BHs comes from studies of X-ray Fe Ka line profiles (Blandford et al. 1990) and by theoretical 
arguments that powerful radio jets are powered by the extraction of energy from rotating 
BHs (Miller 2007). 

In a recent study, Wang et al. (2006; hereafter WCHM) estimated a high average ra- 
diative efficiency of 30%-35% for quasars at 0.4 < z < 2.1, implying that most supermassive 
black holes are rapidly rotating. However, the existence of rotating BHs with 77 > 0.18 is 
not confirmed by magnetohydronamic (MHD) simulations: gas loses more angular momen- 
tum prior to accretion in an MHD disk than in a standard thin disk (Gammie & Shapiro 
2004; Shapiro 2005). Observationally, Shankar et al. (2008) present and review evidence 
for 77 < 0.1. In this work we show that the WCHM method is not independent of quasar 
lifetimes. In § 2 we correct and extend the WCHM method for determining radiative effi- 
ciencies, in § 3 we apply it to a sample of Sloan Digital Sky Survey (SDSS) quasars and in 
§ 4 we discuss our conclusions. 



2. The Method: 

We assume that quasar light derives only from accretion of matter onto a black hole, 
neglecting the effects of BH mergers. Suppose that mass propagates through a thin accretion 
disk around an accreting BH at a rate of M acc during some time interval of At. The BH mass 
growth rate is given by M = (1 — i])M acc , where r\ is the radiative efficiency. The quasar 
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radiates at bolometric luminosity L (Marconi et al. 2004) given by: 



LAt = 7]M acc c 2 At = — - — Mc 2 At. (1) 

1—7] 

We can rewrite the above equation as: 

LAt 5e 

V = = = 7 f— ? 2 

LAt + Mc 2 At 6e + 5pc 2 

where 5e = LAt/V com is the change over the time At in the comoving radiative energy 
density and 5p = MAt/V com is the accompanying change in the comoving BH mass density, 
both due solely to this BH's accretion. 

By analogy, an average radiative efficiency can be defined for any sample of quasars at 
redshift z: 

— Ae(z) 

V{Z) = Ae(z) + Ag(z)c 2 (3) 

where Ae(z) and Ag(z) are the estimated changes in the cumulative radiative energy density 
and the cumulative mass density of the sample in the redshift range (z,z + Az). We now 
define these cumulative densities. 

We require the quasar black hole mass function n(M, z), defined such that n(M, z) AM Az 
is the comoving number density of black holes with masses in the range (M, M + AM) in 
the redshift range (z,z + Az). We also require the quasar bolometric luminosity function 
ijj(L, z), where i/j(L, z) AL Az is the comoving number density of black holes with bolometric 
luminosity in the range (L, L + AL) in the redshift range (z, z + Az). 

Over its lifetime t q , a single quasar accretes a mass Mf ina i=M avg t q and radiates an 
energy L avg t q . We need to express those quantities in terms of the observables Mj and Lj, 
which are the black hole's mass and bolometric luminosity at the redshift of observation, z. 
We do so by defining correction factors / and g such that 

% — 1 J — 1 

where the averages are over all iV quasars in the sample. 

First, consider the mass density that contributes to Eq. [3j Over their lifetimes, all black 
holes observed at redshift z will accrete a comoving matter density of 



^n(Mj, z)M fi na i AM-i = ^n{M h z)M i AM i / f . 

Mi Mi 



(5) 
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The cumulative, lifetime amount of matter accreted by all black holes observed at >z and 
above is the comoving lifetime mass density summed over all redshifts >z: 

g{z) = M^fc) = Azk Yl n ( Mi > ZkWiAMi/f (6) 

z k >z z k >z Mi 

which agrees with Eqs. 2 and 6 of WCHM if / = 1. Notice that units of Ag/Az are mass 
per comoving volume per redshift. 

Now consider the radiative energy component of Eq. [31 Over their lifetimes, all black 
holes observed at redshift z will radiate a comoving energy density of 

^2ip(L j} z)L aV gi q ALj = ^2^p(L j ,z)L j igAL j /g, (7) 

where t q is the average quasar lifetime^ The cumulative, lifetime amount of energy radiated 
by all black holes observed at redshift z and above is the comoving lifetime energy density 
summed over all redshifts above zu 

e{z) = Ae ( z k) = Az kJ2^( L i, zjL&ALj/g (8) 

z k >z z k >z Lj 

which differs from Eqs. 3 and 5 of WCHM. Their expression for Ae(z) is a factor of Atk/ Azki q 
times the true value above (assuming g — 1), where Atk is the cosmological time spanned by 
the redshift interval Azk- The units of As/ Az above are energy per comoving volume per 
redshift, whereas in Eq. 5 of WCHM they are energy per comoving volume per (redshift) 2 , 
which is incorrect. 

The first quantity of interest, the change in the cumulative comoving mass density of 
actively accreting black holes over the redshift range (z, z + Az) — Ag(z) — is Az times the 
sum of the masses of all individual quasar black holes in the sample in that z range, divided 
by the comoving volume: 

Ag(z) =AzJ2 < M i, z ) M i AMi//. (9) 

Mi 



1 Time periods without accretion are not counted in t q . While t q can be defined as the sum of all time 
periods during which an average quasar is actively accreting mass, an acceptable observational definition 
might be the time required to accrete, e.g., 95% of the quasar's final mass (see Hopkins et al. 2006). 

2 For convenience, in Equations 6 and 8 the redshift sum appears in front of the mass or luminosity sum 
for that redshift bin. That latter sum is computed in the feth redshift bin, multiplied by Az k and then added 
to the mass or luminosity sum from the (/c+l)th redshift bin times Az k +i, and so on. The redshift bin size 
does not matter as long as n(M, z k ) or ip(Lj, z k ) does not change considerably within a bin. For example, if 
the redshift bin width was halved, each term in the redshift sum would be half as large but there would be 
twice as many terms, yielding the same result. 
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In other words, the change in the cumulative comoving mass density in the redshift bin 
(z,z + Az), which is Ag(z)/ Az, is the same as the comoving mass density in that bin. 
Similarly, the second quantity of interest, the change over the redshift range (z, z + Az) in 
the cumulative radiative energy density ever observed from all quasars — Ae(z) — is 

Ae(z) = Az ^2^p(Lj,z) Lj ALj t q /g. (10) 



These expressions for Ag(z) and Ae(z) can be substituted into Eq. [3] to find fj(z) for 
the sample of quasars under consideration. Doing so, and grouping i q , f and g together, we 
obtain: 

^2if)(Lj,z) Lj ALj 

fj(z) = * - 2 (11) 

^i L v z ) L i AL i + W ^n(M i} z) Mi A Mi 

where the Az in the numerator and denominator have cancelled out, making the calculation 
independent of Az even if t q is longer than the cosmic time interval corresponding to Az. § 

Thus, the WCHM method for studying radiative efficiencies is not independent of the 
average quasar lifetime t q . A factor of i q enters because both the mass-energy growth and 
the radiative energy output of a quasar must be summed over its entire lifetime (or over the 
same portion of its lifetime). The mass-energy sum yields the final mass-energy of the black 
hole MtotC 2 , while the radiative energy sum yields L avg i q . However, this method remains 
independent of obscured sources and is also powerful because it can be implemented for any 
sample of quasars regardless of selection effects. It requires only that the changes in the 
cumulative mass-energy and radiative energy densities be computed using the same objects. 
Of course, selection effects will determine if the resulting radiative efficiency is relevant for 
quasars in general. 

Also, WCHM in effect assumed f — g — 1. We adopt more realistic estimates by 
examining Fig. 14 of Springel et al. (2005) and Fig. 1 of Hopkins et al. (2005), from which 



3 For example, consider a population of quasars observed at a redshift z\, and a population observed at 
a redshift Z2 > z\. Those quasars are different objects in different regions of the universe, but the low-z 
population might still be the descendents of the high-z population. (If we could watch both regions of the 
universe over cosmic time, we might see the high-z population evolve into the low-z population.) In that 
case, the masses of the quasars would be counted twice, at different cosmic times, but the radiative output 
would also be counted twice, at those same times. The estimates of the comoving mass and radiative energy 
densities would be systematically off, but in such a way that the radiative efficiency calculation would still 
be accurate. 
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we respectively estimate f — 1 and g ~ 10 for those models. Adopting f=l means that a 
quasar is observed when its black hole has accumulated essentially all of its mass, while g=10 
means that the observed bolometric luminosity of a quasar is around 10 times larger than 
the average luminosity of a quasar throughout its entire life. Our value of g is estimated 
by averaging Lj/L ave over the period of time when quasar shows its final burst of activity 
(between 1.2 to 1.4 Gyr in Fig. 1 of Hopkins et al. 2005). 

If WCHM had assumed / = 1 and g = 10, we estimate that their (incorrect) calculation 
would have yielded r] < 0.1 instead of rj ~ 0.32. However, Eq. 11 shows that even with 
realistic / and g values the average quasar radiative efficiency is dependent on the average 
quasar lifetime. We now explore the implications of this dependency. 
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Fig. 1. — Tracks of putative radiative efficiency rj vs. redshift z for different assumed values 
of the quasar lifetime in Gyr, all assuming / = g = 1. The blue horizontal lines show 
rj for the Extreme-Kerr BH limit (EKERR), the Thorne BH limit (THOR), the magneto- 
hydrodynamic BH limit (MHD) and the Schwarzschild BH limit (SCH) respectively. The 
assumed quasar lifetimes are given next to each track, on the left for tracks chosen to match 
a particular radiative efficiency line at low redshift epochs and on the right for high redshift 
epochs. An Extreme-Kerr BH with t q ~ 0.68 Gyr at z ~ 0.9 would violate the maximum 
possible radiative efficiency at z ~ 2. This inconsistency indicates the existence of a hidden 
parameter; namely, Eddington ratio which is considered in Figure [2j 
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3. Application to SDSS Quasars 

3.1. Estimation of Black Hole Masses and Bolometric Luminosity 

Based on reverberation mapping studies of local active galaxies, an empirical scaling 
relationship has been developed to estimate black hole masses by (e.g.) Kaspi et al. (2000), 
Vestergaard (2002) and McLure et al. (2002). The Sloan Digital Sky Survey (York et al. 
2000) Data Release 3 quasar catalog (Schneider et al. 2005) provides us with a large sample 
of quasars at redshifts 0.7 < z < 2.1 which have their Mg II A2800 emission line redshifted 
into the SDSS spectral range. We have estimated black hole masses for 27728 such quasars 
from the dispersion of the Mg II emission line and the continuum luminosity at A=3000 A 
and for two scenarios. Scenario A: we assume that a black hole mass can be estimated from 
a conventional virial relationship as: M BH /M Q = 30.5 [ALa] - 5 ^^ where L x has units of 10 44 
erg s _1 and CMg km s -1 . Scenario B: we follow Marconi et al. (2008) in considering the effect 
of the radiation pressure of the quasar's radiation on black hole mass estimates, namely that 
it reduces the effective gravity on the broad emission line region, yielding narrower lines 
at a given mass. We adopt the BH mass relationship with this radiation correction to be: 
Mbh/Mq = 5.75Ca[\L x ] - 5 crl lg + £ CT [AL A ] where ( a = 2.4 ± 1.5 and log^ = 6.9 ± 0.5 (Rafiee 
et al. 2008, in prep.). 

Furthermore a Malmquist-like bias has been estimated following § 4 of Shen et al. (2008). 
Our black hole mass estimates have been adjusted downward by a mean bias of 0.4 dex for 
scenario A and 0.11 dex for scenario B, arising from the steepness of the BH mass function 
and the scatter in BH mass estimates. 

We have matched our sample with that of Richards et al. (2006a), a homogeneously 
selected and statistically complete sample of 15343 DR3 quasars with redshifts z < 5 drawn 
from an effective area of 1622 deg 2 . This procedure yields a subsample of 6704 quasars. 
The bolometric luminosities of these quasars have been estimated from = C\\L\ with 
C x = 5.15 for A=3000 A following Shen et al. 2008. 

Comoving volumes and luminosity distances have been calculated using a ACDM cos- 
mology with h = 0.71, Q\ = 0.74 and Q m = 0.26 (Spergel at al. 2007). Corrections have 
been made for the limited areal coverage of the Richards et al. (2006a) sample and for the 
5% incompleteness of the SDSS at 0.7 < z < 2.1. For more details, see Fig. 8 of Richards et 
al. (2006a). 
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3.2. Radiative Efficiency and Redshift Binning 



For thin accretion disks, the only disks we consider in this paper, r\ varies from 0.057 
for a Schwarzschild BH (with a* = Jc/GM 2 = 0, where J is the angular momentum) to 
0.42 for an Extreme-Kerr BH (a* = 1). When one takes into account the effect of radiation 
energy captured by the BH (Thorne 1974), the radiative efficiency reaches a maximum of 
~ 0.30 (a* ~ 0.998); we refer to that Thorne BH. It is more realistic to assume 

an MHD disk wherein magnetic turbulence provides a torque to remove angular momentum 
from the inflowing gas (Shapiro 2005); in that case, the maximum radiative efficiency is 
~ 0.18 (a* ~ 0.938)0 

We divide our sample into twelve redshift bins. In each bin we can compute rj for any 
given value of t q . The results are shown in Figured] as tracks of rj(z) for eight different values 
of t q chosen to match the 7] of a Schwarzschild, MHD, Thorne or Extreme-Kerr black hole 
at either the high or low redshift limit of our sample. Figure Q] is the corrected version of 
Figure 2 of WCHM. Both figures show the evolution of r] with z for a flux-limited quasar 
sample, assuming constant i q . However, this figure does not give a complete picture since 
there is another degree of freedom not being considered; namely, the Eddington ratios of the 
quasars. For example, taken at face value, Figure [1] suggests that a quasar with a lifetime 
of 0.68 Gyr can be powered by a Extreme-Kerr BH at z ~ 0.9 while the same quasar would 
violate the maximum possible radiative efficiency at z ~ 2. Since quasars in our sample have 
larger values of the Eddington ratio at z ~ 2 than at z ~ 0.8, the trend in Figure [1] might 
be explained by the importance of the Eddington ratio rather than the redshift. 



We assume the radiative efficiency may be a function of the Eddington ratio T = 
Lboi/LEdd, where L Edd = 1.26 x 10 38 (M BH /M Q ) erg s _1 . In that case, the changes in the 
cumulative comoving mass and energy densities in the Eddington ratio bin (T, T + AT) are: 



4 We relate r\ to a* assuming co-aligned accretion on to rotating BHs. King, Pringle & Hofmann (2008) 
have pointed out that the effective rj{a*) will be different for randomly-aligned accretion. The conversion 
from r\ to a* will differ for each combination of co- and randomly-aligned accretion, but high r\ will always 
require high a*. 



3.3. Eddington Ratio Binning 




(12) 




(13) 
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where n and ip now depend on T as well as z. The bullet subscript denotes quantities binned 
in T as well as z. 

The average radiative efficiency for a given T and z is 

Ae.(2, T) + A£>.(z, T)c 2 

For the same sub-sample used in section 3.2 and for twelve redshift bins, ten mass bins, and 
ten luminosity bins, the radiative efficiency has been estimated for four different Eddington 
ratio bins (Figure |2J), each containing one quartile of the objects. Figure [2] shows that 
radiative efficiency is not a function of redshift but rather of quasar lifetime and Eddington 
ratio. 



4. Discussion and Summary 

Determinations of quasar lifetimes, Eddington ratios and radiative efficiencies are inter- 
related. Given constraints on (or assumptions about) quasar lifetimes, the WCHM method 
can be used to constrain quasar radiative efficiencies and BH spins. ( Without such con- 
straints, the average quasar rj cannot be estimated by this method.) Conversely, the range 
of radiative efficiencies possible for the full range of BH spins can be used to constrain 
average quasar lifetimes, as long as luminous quasars are not powered by radiatively ineffi- 
cient accretion flows (RIAFs; see, e.g., Blandford & Begelman 1999). For example, for the 
rj = 0.065, T ~ 0.4 model of Shankar et al. (2008), we predict i q ~ 140 million years in our 
scenario A or t q ~ 370 million years in scenario B, which could be used as a further test of 
those models in comparison to others. 

Assuming / = 1 and g = 10 (see the end of § 2), quasar lifetimes can be constrained 
according to the Eddington ratio of the quasar. Lifetimes estimated this way are within a 
factor of a few of literature lifetime estimates. For example, for BHs in the mass range of our 
sample (1O 8 <Mb#<1O 1o M ), a lower limit lifetime of 530 million years can be established 
for black holes with 0.02<T<0.38 (Panel 1 of Figured, scenario A) or around 750 million 
years in scenario B for 0.02<T<0.26. This lower limit corresponds to the Schwarzschild 
case, since a rotating black hole at the same T will require a longer lifetime to build up its 
observed mass. This lower limit lifetime is less than a factor of two lower than the mean 
lifetime of one billion years estimated by Marconi et al. (2004) for T = 0.1 and rj = 0.04 in 
the same range of Mbh- As another example, a luminous quasar powered by a relatively 
low-mass black hole — which may be a typical early stage in a quasar's evolution — will 
have T > 0.4 and can have a typical lifetime of 140 to 510 million years (Panel 4 of Figure 
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[2^i). This range is only a factor of ~ 3 larger than the mean lifetime of 30 — 130 million years 
estimated by Yu & Tremaine (2002) for luminous quasars, and is consistent with the mean 
lifetime of 100 — 450 million years estimated by Marconi et al. (2004) for super-Eddington 
accretors. 

In principle, given constraints on T a r\M j (1 — T])Mbh, f/g and t q for quasar samples, 
one could estimate the historical frequency of RIAF episodes in those quasars by plotting 
quasar lifetimes versus Eddington ratios. For example, consider quasars lying below the 
Schwarzschild curves in Figure [3] (the normalization of which is a function of f/g, as seen 
by comparing the two axis in either scenario A or B). Above our lower mass limit of 10 8 M Q , 
such quasars must either have had a RIAF phase in order to explain their observed masses, 
or they must have observed Eddington ratios lower than their historical average: Y f, s < T. 
(In the latter case, the quasars historically would have been located horizontally to the right 
in the diagram, lying between the Schwarzschild and Thorne curves at a value of T = T 
sufficient to yield the observed Mbh in the observed i q .) Conversely, quasars lying above the 
Thorne curves in Figure [3] require T Q b s > T. A low T might result if the BH spin does not 
increase as fast as the BH mass does, perhaps due to counter-rotating gas accretion phases. 
If / = 1 and g = 10 and t q < 1 Gyr, and if RIAFs are unimportant, then not many BHs 
with T < 0.2 can be rapidly spinning. On the other hand, if / = 1 and g = 10 and t q < 140 
Myr, then RIAFs or BH mergers must be important for quasars regardless of T, since only 
then could the observed masses be reached in the inferred lifetimes. 

What can we conclude if we assume that the Marconi et al. (2008) correction for the 
effect of radiation pressure on quasar BH masses is valid, and that t q ^ 1 Gyr (Martini et 
al. 2004; Marconi et al. 2004)? First of all, most quasars can not be rapidly spinning in 
that case. Second, most of the quasars in our sample have a radiative efficiency of < 0.14 
consistent with the results of Yu & Tremaine (2002). This 77 being lower than the MHD 
prediction of Shapiro (2005) might be explained by the effects of BH mergers or by the 
fraction of maximally spinning BHs being low, at least in our sample. 

Alternatively, if one assumes thick disk accretion, where the relation between the BH 
spin and radiative efficiency differs from thin disk accretion, then radiatively inefficient ac- 
cretion becomes more important even for MHD accretion and spinning BHs, making a lower 
i] more plausible. 

In conclusion, the Wang et al. (2006) method, despite its advantages, can only estimate 
the radiative efficiency of quasars and ultimately the spin of black holes if we know enough 
about the accretion process and the evolutionary history of black holes. Better estimates 
of /, g and i q from ultimately, however, more detailed evolutionary models, might improve 
the reliability of the results. Lack of knowledge about the geometry and dynamics of the 
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accretion disk limits the level of reliability of this method. 
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Fig. 2. — Radiative efficiencies vs. z for different Eddington ratio bins and constant assumed 
quasar lifetimes in Gyr (printed in each panel for each BH type), (a) Scenario A. (b) Scenario 
B. Lifetimes were chosen so the resulting r\ would match that of each BH type. We assume 
/ = 1 and g = 10 — quasars are observed at their final mass and ten times their average 
luminosity The shaded areas are the regions corresponding to i q < 1 Gyr, the possible upper 
limit for quasar lifetime suggested by Marconi et al. (2004). 
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Fig. 3. — Quasar lifetime versus Eddington ratio, (a) Scenario A. (b) Scenario B. Quasars 
will lie in the area between the Schwarzschild and Thorne curves if they are described by 
radiatively efficient accretion disks. Below the Schwarzschild curves requires radiatively 
inefficient accretion or T > T & s , while above the MHD curves requires T < T t> s . Left axis 
labels are for f = g = 1; right axis labels are for g = 10 and / = 1. 



